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What is the quantity and composition of material in the Universe? This is one 
of the most fundamental questions we can ask about the Universe, and its answer 
bears on a number of important issues including the formation of structure in the 
Universe, and the ultimate fate and the earliest history of the Universe. Moreover, 
answering this question could lead to the discovery of new particles, as well as 
shedding light on the nature of the fundamental interactions. At present, only a 
partial answer is at hand: Most of the material in the Universe does not give off 
detectable radiation, i.e., is “dark;” the dark matter associated with bright galaxies 
contributes somewhere between 10% and 30% of the critical density (by comparison 
luminous matter contributes less than 1%); baryonic matter contributes between 
1.1% and 12% of critical. The case for the spatially-flat, Einstein-de Sitter model 
is supported by three compelling theoretical arguments — structure formation, the 
temporal Copemican principle, and inflation — and by some observational data. 
If Q is indeed unity — or even just significantly greater than 0.1 — then there is a 
strong case for a Universe comprised of nonbaryonic matter. There are three well 
motivated particle dark-matter candidates: an axion of mass 10 -6 eV to 10 -4 eV; 
a neutralino of mass 10 GeV to about 3TeV; or a neutrino of mass 20 eV to 90 eV. 
All three possibilities can be tested by experiments that are either being planned 
or are underway. 


I. Weighing the Universe: Dark Matter Dominates! 

The Friedmann-Robertson-Walker cosmology, also known as the hot big bang model, 
provides a reliable and tested accounting of the Universe from about 10 -2 sec after the 
bang until the present. It is so successful that it is known as the standard cosmology. In 
the context of this cosmology the critical density separates models that expand forever 
\ ( p < pcRlT ) from those that ultimately recollapse ( p > pcrit); Pcrit = 3Hq/8ttG cz 
i 1.8S h I. 2 * x 10 -29 gcm -3 ~ 1.05/i 2 x 10 4 eVcm -3 , where the present value of the Hubble 

■ parameter Ho = 100/ikmsec -1 Mpc -1 ~ 1 / 3000A - 1 Mpc. I will denote the ratio of the 
1 total energy density p (including a possible vacuum energy) to the critical density by Cl = 
1 p/pcRIT, and the fraction of critical density contributed by species i by, fl, = pi/pcRlT- 
j The flat Einstein-de Sitter model corresponds to fi = 1; the negatively curved model to 

■ Q < 1; and the positively curved model to Q > 1. The radius of curvature can be expressed 


in terms of H 0 and ft: Rcurv = H$ x /\£l — l) 1 / 2 . 

There are a variety of methods for determining ft . 1 Broadly speaking they can be 
divided into two qualitatively different categories. First, there are the dynamical methods 
where the mass density is inferred by its gravitational effects; these include measuring 
the “rotation curves” of spiral galaxies, the virial masses of clusters of galaxies, and the 
local peculiar-velocity field. Second, there are the kinematic methods, which are sensitive 
to both the space-time geometry and the time evolution of the cosmic scale factor R(t). 
They include the classic Hubble diagram (red shift-luminosity relation), the red shift- 
galaxy count relation, red shift-angular size relation, and others. 2 

Dynamical Methods 

One can use Kepler’s third law to determine the mass of a galaxy: GM = v 2 r, where v 
is the orbital velocity of a “test particle,” r is its orbital radius, and M is the mass interior 
to the orbit (valid for a spherical mass distribution); or its statistical analogue, the virial 
theorem, to determine the mass of a gravitationally bound cluster: GM = (v 2 )r where 
M is the cluster mass, { v 2 ) 1 / 2 is the velocity dispersion of the galaxies, and r is the core 
radius of the cluster (orbits are assumed to be distributed isotropically). 

For simplicity, one can imagine that one uses these methods to determine the “average 
mass per galaxy” and then multiplies it by the number density of galaxies to determine 
the average mass density p. In reality, astronomers use these methods to determine the 
mass-to-light ratio for spiral galaxies and for clusters of galaxies; from the mass-to-light 
ratio they infer the average mass density 

p=[M/L)C, (1) 

where (M/L) is the mass-to-light ratio, and C is the luminosity density, whose value 
is about 2.4/i x 10 8 TsgMpc -5 in the Bt system. The critical mass-to-light ratio is 
(M/Z-)crit — 1200 h Mq/ Lq, where subscript 0 refers to solar units. 

“Rotation curves” — that is orbital velocity as a function of orbital distance — have been 
determined for numerous spiral galaxies. They are obtained by measuring the Doppler 
shifts of stellar spectral features and of the 21 cm radiation from neutral gas clouds (HI 
regions) — the stars and clouds act as gravitational test particles. Rotation curves axe 
all qualitatively similar; they rise rapidly from the galactic center and remain flat (v — 
const) out to the furthest distances that can be probed — eventually, one “runs out” of 
test particles, i.e., stars and gas clouds. Since v = const implies M(r) oc r, this means 
that one “runs out” of stellar light and 21 cm radiation before the mass of the galaxy has 
“converged.” In some cases the 21 cm rotation curves have been determined to a distance 
that is three times that where the light has fallen to 1% of its value at the center of the 
galaxy. 

By restricting oneself to the bright central regions of a galaxy one can use the rotation 
velocity to infer the amount of mass associated with the “luminous” part of the galaxy; 
doing so one finds that luminous matter contributes 

ftLUM ~ 0-01, (2) 

which is far from the critical density. A similarly small value is obtained by using the 
mass-to-light ratio determined for the local solar neighborhood, (M/L) i OC al ~ 2 — 3. 
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Fig. 1 . Upper: F-band surface brightness of NGC 3198 in units mag axes ec 2 (F-band 
covers a “red” part of the spectrum from about SOOOA to about 7000A). 21 cm rotation 
curve for NGC 3198 (dots with error flags) and rotation curve predicted from the luminous 
matter alone (assuming constant mass-to-light ratio M/Lq = 4). Lower: Rotation curves 
for a number of spiral galaxies determined from 21 cm observations. Vertical bars indicate 
the point beyond which the surface brightness is less than 25 (blue) mag arcsec“ 2 (less than 
about 1% of the central surface brightness). [From Sancisi and van Albada in Kormendy 
and Knapp, Ref. 1.] 





Based upon the fact that many rotation curves stay flat out to distances far beyond 
where the surface luminosity of the galaxy is negligible, one can infer that there is much 
more matter associated with spiral galaxies that is dark (i.e., does not give off visible 
radiation) than is luminous. For our own galaxy the rotation velocity has been measured 
out to a distance of about 20 kpc, at which point the dark matter contributes about three 
times more mass than the luminous matter (for reference the solar system is about 8.5 kpc 
from the center of the galaxy). There is weaker evidence that this dark matter exists in a 
spherically-symmetric, extended halo with a density that varies as r -2 at large distances 
from the center of the galaxy. 

Based upon the rotation curves, one can conclude that the dark halo material in spiral 
galaxies contributes at least three to ten times the mass density that luminous matter does, 

^HALO >0.03-0.10. (3) 

Since there is no convincing evidence for a rotation curve that “turns over” and decreases 
as r -1 / 2 indicating that the halo mass has converged, it is possible that the halos of spiral 
galaxies extend a factor of order ten further and thereby provide the critical density." 

The existence of dark matter halos in spiral galaxies provides the answer to one puzzle — 
the stability of galactic disks — and raises another — the apparent conspiracy of the luminous 
matter and dark matter to produce smooth rotation curves. A disk-like structure is subject 
to many instabilities, and a massive halo stabilizes a disk-like structure against these 
instabilities thereby resolving a longstanding puzzle. However, the existence of dark matter 
halos raises another question: Why do the inner and outer parts of the rotation curve join 
so smoothly, in light of the fact that the inner part of the rotation curve is supported by 
luminous matter and the outer part by dark matter? (The rotation curves of most spiral 
galaxies are very similar, with the rotation velocity rising rapidly from zero at the center 
to a nearly constant value. The rotation curve for our own galaxy is quite flat and smooth 
at our position, in spite of the fact that the gravitational support for rotation velocities 
at our position are about equally split between luminous disk material and dark halo 
material.) Some (e.g., Peebles) have argued that this is evidence that the halo and disk 
have a similar composition — baryons, while others put their faith in numerical simulations 
of the formation of galactic halos and disks that indicate that this occurs quite naturally 
when the ratio of nondissipative dark matter and dissipative luminous matter is of order 
ten. 

There is some evidence that individual elliptical galaxies contain significant amounts 
of dark matter, although the case is not as well established as that for spirals. Most cluster 
galaxies are ellipticals, and as I will now discuss there is strong evidence for dark matter 
in clusters. 

Estimates of the mass density based upon the virial masses of clusters lead to 

^CLUSTER — 0.1 - 0.3, (4) 

° There are arguments to the contrary; e.g., mass estimates of the Milky Way and 
Andromeda based upon their velocity of approach seem to indicate that their halos could 
not be this large, although such arguments assume that the Milky Way and Andromeda 
are on a radial orbit and are approaching each other for the first time. Likewise, mass 
estimates of the Milky Way based upon the orbits of its satellite galaxies indicate the 
same, although it is assumed that the orbits are isotropically distributed. 3 
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again indicating substantially more mass than that required to account for the light. Sev- 
eral points should be noted: (1) X-ray emission from hot intracluster gas indicates the 
presence of comparable or greater amounts of baryonic mass than that associated with 
the visible light (dark is a relative term!), but no where near enough to account for the 
cluster’s virial mass. (2) Since only about one in ten galaxies resides in a large cluster, one 
can question whether or not the mass-to-light ratio — and value of 0 — deduced from clus- 
ters is indicative. However, there seems to be no question that clusters contain significant 
amounts of dark matter. (3) These determinations are based upon the assumption that 
the clusters axe well virialized, single objects and that the galaxy orbits are distributed 
isotropically; moreover, the cluster core radius is inferred from the distribution of the 
visible galaxies. If galaxies have sunk deep into the cluster potential, e.g., due to dynami- 
cal friction, then the actual core radius of the cluster — and cluster mass — could be much 
larger 30 (just as with galactic halos). 

Peculiar Velocities 

The velocity of a galaxy can be split into two pieces: the velocity due to the general 
expansion of the Universe (or Hubble velocity) which is radial and proportional to galaxy’s 
distance from us; and the peculiar velocity, the velocity the galaxy has in addition to its 
Hubble velocity. 0 Any peculiar velocity that is not “supported” by a gravitational field will 
decay with time, inversely with the cosmic scale factor R(t). Put another way, peculiar 
velocities arise due to the lumpy distribution of matter — and thereby offer a probe of the 
density field. In contrast, the distribution of bright galaxies only probes the distribution 
of light — and the two distributions need not be the same. 

In the linear perturbation regime, i.e., Sp/p < 1, the Fourier expansion of the velocity 
field, Vfc, is related to that of the density field, 8k, v* = — tkJ?(t)<jfc(t)/|k| 2 , and to a good 
approximation |v*| ~ fl° 6 Ho\Sk\/k. Suppose the peculiar velocity of an object is primarily 
due to linear perturbations on the scale A, then 



Even if the contribution from one Fourier component does not dominate, Eq. (5) still 
illustrates the correct dependence of the peculiar velocity upon Cl. d One can exploit this 
relationship in different ways: (i) input and Sv to infer Sp(r)/p; (ii) input and 6p(r)/p 
to infer 6v; or (iii) input <5v and 6p(r)/p to infer Cl. The last of these alternatives is the 
one we are interested in here; however, what one can directly measure is ^(rj/no, and 
so one must relate Sno/n g to Sp/p (uq is the number density of “bright” galaxies). The 
simplest ansatz is to take them to be equal: “light traces mass.” A slightly more general 

b It should be mentioned that in 1933 the astrophysicist Fritz Zwicky pointed out that 
the mass associated with the light in several clusters was much less than the mass required 
to bind the cluster — and thus was the first to identify the dark matter problem. 

c Of course, we can only measure the component of the peculiar velocity that is parallel 
to the line of sight. 

d More precisely, the peculiar velocity at position r is £v(r) = —CI°- 6 (Ho/4tt) f 6p( r')(r — 
r')dV/|r - r'| 3 p. 
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approach is to assume that “light is a biased tracer of mass:” Sp/p = b~ 1 (Snc/nG), where 
1 < b < 3 is the biasing factor. 

Using the IRAS catalogue of infrared- selected galaxies to determine the mass distribu- 
tion (i.e., Sno/nc), several groups have used measurements of the local peculiar-velocity 
field 4 to infer f l 0 6 /b ~ 1., with an estimated uncertainty of about 0.3 or so. 5 With some 
delight, I note that this technique seems to suggest that ft is indeed close to unity. Al- 
though I caution the reader that these results are still preliminary, if they hold up, they 
will provide the strongest evidence to date for a large value of ft!* 

Before going on to the kinematic methods, I mention that there are other dynamical 
methods for determining ft, including the use of gravitational lens systems to measure 
cluster and galaxy masses, Virgo infall (which is similar to the peculiar-velocity method 
mentioned above), cosmic virial theorems, and pair-wise velocities of galaxies. 6 

In addition, there may or may not be another, more local dark matter problem. The 
mass density of the disk in our neighborhood can be determined by studying the motions 
of stars perpendicular to the plane of the disk, and by a “direct inventory” of the material 
in the local neighborhood (stars, white dwarfs, gas, dust, etc.). In principle the two results 
should agree. The local mass density inferred from dynamics/ 1.3 x 10 -23 g cm -3 , is about 
a factor of two larger than can be accounted for by the local inventory. 7 This discrepancy 
of a factor of two may or may not be significant. In any case, it has little bearing on the 
“big” dark matter problem. Since the mass density of the local neighborhood is dominated 
by luminous matter, this additional dark matter — If it exists — makes a contribution to ft 
that is at most comparable to that of luminous matter. 

Moreover, this local dark matter cannot be due to halo material: Based upon the 
rotation curve of our galaxy and detailed models for the distribution of matter in our 
galaxy, the local halo density is estimated to be 8 

Phalo — 5 x 10 -25 g cm -3 ~ 0.3 GeV cm -3 , (6) 

with an uncertainty of about a factor of two. The local halo density is about a factor of 
ten smaller than the local disk dark-matter density; put another way, if the halo material 
accounted for the disk dark-matter density, the local rotation velocity would be about a 
factor of three larger than its measured value! 

Kinematic Determinations 

There are a number of classic kinematic tests — luminosity-red shift (or Hubble dia- 
gram), angle-red shift, galaxy-number count -red shift — that can in principle be used to 
determine our cosmological model. 2 These tests depend upon the global space-time geom- 
etry and the time evolution of the scale factor. For example, the luminosity distance to 

* The infrared bright galaxies tend to be spiral galaxies in the field, and so clusters are 
under-represented. The authors have tried to correct for this by including some important 
clusters, and find that their results do change significantly. One must also worry about 
convergence; that is, has one reached the point where the contribution of galaxies at still 
larger distances has become insignificant. 

f This density is known as the Oort limit, in honor of the first astronomer to address 
this problem. 
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a galaxy at red shift z, d\ = C/AirJ-, is related to the coordinate distance to the galaxy, 
r(z), by 

d L = r ( 2 ) 2 (! + *) 2 , 



where the present value of the scale factor R(to) is taken to be one and (1 + 2 ) = R(t)~ l . 
Since the evolution of the scale factor depends the equation of state, e.g., p = 0, matter- 
dominated, R oc t 2 / 3 ; p = p/ 3, radiation-dominated, R oc t 1 / 2 ; p — —p, vacuum- 
dominated, R oc exp (Ht), the functional dependence of r(z ) does too. Thus, the red 
shift -luminosity distance relation depends upon both the curvature of space and the com- 
position of the Universe. For a matter-dominated model 

Hq(1l = q 0 2 zqo + (<70 ~ 1 ) + 1 — l^j = z [1 + (1 — 50 )^/2 + ••■]■ ( 8 ) 

where qo = —Rq/Hq = ft(l -f 3p/p)/2 = ft/2, and the second expression is an expansion 
in 2 . 

The success or failure of this technique depends upon obtaining accurate luminosity 
distances for objects out to red shifts of order unity. Accurate luminosity distances requires 
the existence of objects of known luminosity (standard candles). Here lies the problem; 
evolutionary effects are likely to be important, especially at high red shifts, and it is 
difficult to determine even the sign of the evolutionary effects let alone reliably estimate 
the magnitude! Nevertheless, there axe some who believe that the I \ -band ( 2.2 pm) version 
of the Hubble diagram will prove useful , 9 as evolutionary effects are lessened . 5 

A kinematic test with great cosmological leverage and promise is the galaxy count -red 
shift relation. The number of galaxies seen in the red shift interval dz and solid angle 
du depends upo n the nu mber density of galaxies na(z) and the spatial volume element, 
d,V = r 2 drdu /\/ 1 — fcr 2 . This relationship too depends upon both the spatial curvature 
and the time evolution of the scale factor. For a matter-dominated model, 

dA^GAL _ ^GAl ( z )[ z ?0 + (go ~ 1 )(v/ 2 < 7 qZ + 1 — l)] 2 
dudz dfg(l + z) 3 9 q [1 — 2go + 2go(l + 2)] 1 / 2 

— ^ 2 ^gal(- ? )[ 1 — 2(g 0 + 1)2 H ]/H$. (9) 

For fixed (comoving) number density of galaxies, the galaxy count increases with decreasing 
(or qo) because of the increase in spatial volume. Loh and Spillar 10 have used the galaxy 
count-red shift test with a sample of about 1000 field galaxies — red shifts out to 0.75 — 
to infer ft = O. 9 I 0 5 (95% confidence). Their result has drawn much criticism; in part 
because their red shifts are not spectroscopically determined (they are determined by six- 
band photometry) and because their results are sensitive to the assumptions made about 
galactic evolution . 11 

9 When one observes a galaxy of moderate red shift in the visible, the light one sees 

comes from the blue or UV part of the spectrum and is produced by massive stars that 
evolve rapidly. By contrast, observing in A'-band, the light one sees was emitted in the red 

part of the spectrum and is produced by lower mass stars that evolve much more slowly. 
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Fig. 2. Predicted light -element abundances as a function of the baryon-to-photon ratio 
■q in the standard scenario of big-bang nucleosynthesis; error flag indicates the change 
in 4 He that arises for Arj^(n) — ±0.2 min. The inferred primordial abundances are: 
Yp = 0.24 ± 0.01; D/H > 10" 5 ; (D± 3 He)/H < 1.1 x 10" 4 ; and 7 Li/H ~ 1.2 ± 0.3 x 10" 10 . 
Concordance between the predicted and measured abundances requires: 3 x 10“ 10 < 77 < 
5 x 10" 10 ; or 0.011 < 0.011ft" 2 < Q B < 0.019ft" 2 < 0.12. 




In principle, this test is less sensitive to evolution, provided that the number of galaxies 
remains constant and their luminosities do not evolve so drastically that they cannot be 
seen. Recent deep galaxies counts indicate an excess of galaxies at higher red shifts — 
indicative of a low value of fi. 12 (If galaxy mergers axe very important — as they may well 
be in cold dark matter scenarios — the number density of galaxies at higher red shifts would 
be expected to be larger.) At the moment, determinations of 0, based upon the galaxy 
number count test axe not conclusive. However, many believe that this method has great 
potential because a large sample of objects can be used and it is less sensitive to evolution. 

Primordial Nucleosynthesis and fie 

Primordial nucleosynthesis provides the most stringent and earliest test of the stan- 
dard cosmology, probing it back to the epoch when T ~ MeV and t ~ sec. The pri- 
mordial abundances of D, 3 He, 4 He, and 7 Li predicted in the standard (and simplest) 
model of primordial nucleosynthesis agree with the inferred primordial abundances of 
these light elements. 13 Moreover, this agreement can be used to constrain one cosmological 
parameter — the baryon-to-photon ratio rj — and one parameter of the standard model — the 
number of light neutrino species N v . 14 Concordance between theory and observation re- 
quires: 

3 x 1(T 10 < tj < 5 x HT 10 and N v < 3.4 

The constraint to the number of light neutrino species has recently been confirmed by 
precise measurements of the properties of the Z° boson, 15 which imply N„ = 3.0 ± 0.1. 
This is an impressive confirmation of the standard cosmology at this very eaxly epoch. 

Primordial nucleosynthesis provides the most precise determination of the baryon den- 
sity. In converting the baryon-to-photon ratio to the fraction of critical density contributed 
by baxyons two other parameters are needed: (i) the temperature of the cosmic microwave 
background (CMBR), which is now accurately determined to be 2.736 ± 0.01 K; 16 and 
(ii) the not so well known value of the present Hubble parameter, 0.4 < h < 1.0 . h The 
nucleosynthesis constraint can be written as 

0.011 < 0.011/r 2 < n B < o.oi9r 2 < 0.12. (10) 

Summary of Our Knowledge of Q. 

What then is the present state of our knowledge concerning the mass density of the 
Universe? Let me try to summarize: 

• Luminous matter contributes only a small fraction of the critical density: Qlum ~ 0.01. 

• Based upon primordial nucleosynthesis baryonic matter contributes: 0.011 < < 

0.12.* 

• Based upon dynamical methods, the mass density associated with bright galaxies is 
f^ABG — 0.2 ±0.1 (the ±0.1 is not meant to be a formal uncertainty estimate). 

h It also assumed that the only change in the baxyon-to-photon ratio since the start of 
nucleosynthesis is the factor of 4/11 decrease caused by the transfer of the entropy in e ± 
pairs to photons when T ~ m e /3. 

' If is close to unity and the cosmological constant is zero, then h must be close to 0.5 
to insure a sufficiently elderly Universe; in this case: 0.04 < Qq ^ 0.12. 
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• There is some evidence that fl might be close to unity; e.g., analyses of the local 
peculiar-velocity field based upon the IRAS catalogue of galaxies, and the result of 
Loh and S pillar. 
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From this I would make the following inferences: 

★ The dark component of the mass density dominates the luminous component by at 
least a factor of ten, and closer to a factor of 100 if = 1, and is more diffuse than 
the l umi nous component, e.g., the halos of spiral galaxies. 

* There is strong evidence for the existence of a dark component of baryons. This should 
not be too surprising since baryons can exist in a variety of low luminosity objects — 
white dwarfs, neutron stars, black holes, brown dwarfs, jupiters, etc. 

★ At present there is no irrefutable case for a universal mass density that is larger than 
that permitted for baryons. 

* If Q is significantly greater than 0.1 — which is already suggested by mass-to- light ratios 
determined for clusters and the local peculiar velocity field — then there is a strong case 
for nonbaryonic dark matter. As I will discuss, there are three attractive particle dark- 
matter candidates whose relic abundance is expected to be close to critical: the axion, 
the neutralino, and a light neutrino. 

* If fl is one, a discrepancy must be explained: why the estimates for the amount of 
material associated with bright galaxies is a factor of about five smaller. There are 
two possibilities. The first, as previously mentioned, the halos of spiral galaxies could 
extend far enough to account for Q, = 1 (and likewise for clusters). Second, there could 
be a component of the mass density that is more smoothly distributed, contributes 
(Ism — 0.8, and is not associated with bright galaxies; e.g., a population of low- 
luminosity galaxies that is more smoothly distributed than the bright galaxies — so- 
called biased galaxy formation — or a relic cosmological constant (more later). 

★ There may be several dark matter problems — and with different solutions. While the 
most economical approach is to assume that all dark matter has the same composition, 
that need not be the case. As mentioned above there is already evidence that some of 
the baryonic matter is dark. Moreover, if there is indeed a local dark matter problem, 
its solution must involve “particles” that can dissipate energy and condense into the 
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disk; it is very unlikely that axions, neutralinos, or neutrinos can do so. Taken at face 
value the observations seem to indicate that there is more dark matter in clusters (per 
galaxy) than in the halos of spiral galaxies — and if ft — 1 — even more dark matter that 
is not associated with clusters. 

To give a concrete example, consider an ft = 1, neutrino-dominated Universe (m v ~ 
92 h? eV). Because of their high speeds, neutrinos would be unlikely to find their way 
into potential wells as shallow as those of galaxies or perhaps even clusters. They 
would likely remain smooth on scales up to the neutrino free-streaming length, Aps — 
40Mpc/(rn„/30eV). The dark matter in galaxies would be baryons — perhaps white 
dwarfs that formed relatively recently in the local neighborhood and brown dwarfs that 
formed when the galaxy did in the halo — and the dark matter in clusters would be the 
neutrinos that eventually made their way into clusters. J 

A Theoretical Prejudice 

While the hard observational evidence for the flat, Einstein-de Sitter model is less 
than overwhelming, there are several compelling theoretical arguments: (i) the temporal 
Copernican principle — if ft ^ 1 the deviation of ft from unity grows as a power of the scale 
factor, begging one to ask why ft is just now beginning to differ from unity; (ii) structure 
formation — in ft < 1 models there is less time for the growth of density perturbations and 
larger initial perturbations are required; in fact, ft < 0.3 models with adiabatic density 
perturbations are inconsistent with the isotropy of the CMBR (see Bond’s contribution to 
these proceedings); and (iii) the flat, Einstein-de Sitter model is an inescapable prediction 
of inflation. To be sure, these arguments are not rooted in hard facts; however, the are 
sufficiently compelling to create a strong theoretical prejudice for ft = 1. From this point 
forward I will adopt this prejudice! 

Dark Matter: New Physics or New Particles 

Finally, there are some who have suggested another explanation for the dark matter 
problem: A deviation from Newtonian (Einsteinian) gravity at large distances. 18 Newto- 
nian gravity (i.e. , the weak field, slow velocity limit of general relativity) is well tested at 
distances from order 10 2 cm to the size of the solar system, order 10 14 cm. However, the 
dark matter problem involves distance scales of order 10 23 cm and greater. If gravity were 
for some reason stronger on these scales there would perhaps be no need for additional 
“unseen” matter to explain flat rotation curves. For example, if G were a function of 
distance, say G(r ) oc r, then flat rotation curves would be consistent with constant mass 
interior to r — eliminating the need for unseen matter. 

I opt for unseen matter. First, it seems unlikely that the same functional dependence 
for the strength of gravity could fit all the observations: While all spiral galaxies have flat 
rotation curves, the size of the luminous part of the galaxy can vary by almost a factor of 
ten, and clusters are even larger. Perhaps a more important reason is that of aesthetics: 
Not only is there no theoretical motivation for such a theory, but it seems difficult, if not 
impossible, to construct a relativistic theory of gravity in which G increases with distance. 

} In a neutrino-dominated Universe it is probably necessary for the dark matter in galax- 
ies to be baryonic, as there seems to be evidence for dark matter in several dwarf galaxies 
in which there is not enough phase space to contain the necessary numbers of neutrinos. 17 
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The one such theory I am aware is extremely complicated and leads to an unsatisfactory 
cosmology. 18 Were it the other way around — lack of compelling dark matter candidates and 
an attractive alternative theory of gravity — I would opt for new physics in the gravitational 
sector. 

II. Why Not Baryons? 

Given the existing observational evidence one has to be bold to insist that fi = 1. 
Moreover, this assumption seems to require one to go still further and postulate that most 
of the matter in the Universe is comprised of particles whose existence is still hypothetical! 
Before taking the big leap, I will comment on the possibility that baryons could contribute 
the critical density. There are two obstacles to this possibility: the nucleosynthesis con- 
straint, ~ 0.12; and finding a place to hide the more than 99 invisible baryons for every 
visible baryon. 

A number of different schemes have been suggested to evade the nucleosynthesis bound, 
for example, massive relic particles that decay into hadrons shortly after nucleosynthesis 
and initiate a second epoch of nucleosynthesis. 19 This scenario requires an unstable particle 
species with very special properties, and seems to lead to the overproduction of 6 Li and the 
underproduction of 7 Li. Perhaps the most clever idea is the scenario where the baryon-to- 
photon ratio is reduced after nucleosynthesis because photons suddenly come into thermal 
contact with “shadow particles” at a lower temperature, which leads to entropy transfer 
from the photons to the shadow world. 20 

Inhomogeneous Nucleosynthesis 

The alternative to the standard scenario that has attracted the most attention is in- 
homogeneous nucleosynthesis. 21 If the quark/hadron transition is strongly first order and 
occurs at a relatively low temperature (< 125 MeV), baryon number can become concen- 
trated in regions where the quark-gluon plasma persisted the longest. Moreover, due to the 
difference in the mean free paths of the proton and neutron around the time of nucleosyn- 
thesis, the high baryon density regions will become proton rich. Clearly, nucleosynthesis 
proceeds very differently, and two new parameters arise: the density contrast between the 
high and low baryon density regions and the separation of the high density regions. 

While early calculations, done with two independent “zones” of differing baryon num- 
ber density and proton fraction, suggested that Qb ~ 1 could be made consistent with the 
observed light element abundances by an appropriate choice of these two parameters, more 
detailed calculations that allow for diffusion between the zones indicate that the predicted 
abundances for all four light elements conflict with observations if Hq ~ 1 — for all values 
of the two parameters. 22 While this appears to be a sad end to an interesting idea, it does 
serve to emphasize the brilliant success of standard nucleosynthesis: The simplest model 
with no extra dials or knobs correctly predicts the primordial abundances of D, 3 He, 4 He, 
and 7 Li. 

Where Is It ? 

Should one be able to evade the nucleosynthesis bound the next problem that one faces 
is where to put all those dark baryons. Ordinary stars, dust, and gas would all be “visible” 
in one way or another. Black holes and neutron stars do not necessarily provide an easy 
way out either. If, as seems likely, black holes and neutron stars evolve from massive 
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stars, where axe the heavy elements these stars produced? And remember, one is trying to 
hide 99 baryons for every baryon that is in a star. Perhaps massive black holes can form 
without overproducing heavy elements; however, there are other worries. If these black 
holes are too massive they will puff up the disk of the galaxy and disrupt binary stars by 
their gravitational effects, and lead to the (unobserved) lensing of distant QSOs. These 
considerations restrict the mass of black holes in the halo to be less than about 10 5 Mq. 23 

White dwarfs, brown dwarfs (stars less massive them about 0.08M© which do not get 
hot enough for to burn hydrogen), or jupiters are better candidates. 24 All could have 
escaped detection thus far and might be detectable in planned experiments to look for 
microlensing of stars in the LMC by such objects in the halo of our galaxy. However, there 
is the issue of the large number of these objects needed. When one smoothly extrapolates 
the observed IMF (initial mass function of the most recent generation of stars) to these 
very small masses, one concludes that are far too few of these objects to account for the 
dark matter in the halo. It should be noted that the IMF is an empirical, rather than 
fundamental, relation, and some have suggested that when the galaxy formed most of its 
mass could have fragmented into small objects. 

To summarize, it is not impossible to evade the nucleosynthesis bound, and there 
is no devastating argument to preclude astrophysical objects comprised of baryons from 
contributing critical density. However, the elegance of the nucleosynthesis argument and 
the difficulty of hiding so many baryons seem to suggest that nonbaryonic dark matter is 
a more promising option to pursue! 

III. Particle Dark Matter 

According to the standard cosmology, 25 at times earlier than the epoch of matter- 
radiation equality, t < *eq = 4.4 x 10 10 (Qh 2 ) -2 sec and T > Teq = 5.5 (Q,h 2 ) eV, the 
energy density of the Universe was dominat ed by a ther mal bath of particles at temperature 
T. For reference, for t < q, T ~ GeV/ y/t / 10 6 sec. 

While the extrapolation of the standard cosmology to very early times (f « 1 sec) 
is a bold step, there are several reasons to expect that such an extrapolation is at least 
self consistent, if not correct: (1) The splendid success of big bang nucleosynthesis, which 
tests the standard cosmology well into its radiation-dominated phase; (2) The fact that 
according to the standard model of particle physics the fundamental degrees of freedom are 
pointlike quarks and leptons, gauge bosons, and Higgs (scalar) bosons* whose interactions 
are expected to remain perturbatively weak at very high energies; and (3) Quantum correc- 
tions to general relativity should be very small for times t 10 -43 sec and temperatures 
T <C 10 19 GeV. 

The implications of this hot, early epoch for cosmology, and dark matter in particular, 
are manifold: At temperature T all particles of mass less than T should be present in num- 
bers comparable to that of the photons; several phase transitions should take place (e.g., 
quark/hadron transition, chiral symmetry restoration, and electroweak symmetry restora- 
tion); and in the symmetry restored phase, the strength of all interactions — including “very 
weak” interactions that have yet to be discovered — should be comparable. 

* Of course, the existence of the Higgs sector has yet to be confirmed, and there could 
well be some surprises at energies greater than 1/y/Up ~ 300 GeV, corresponding to times 
earlier than 10 -11 sec. 
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While the standard SU(3)c ® SU(2)i ® f/(l)y gauge theory of the strong and elec- 
troweak interactions does not offer any dark matter candidates — beyond the now dim hope 
that inhomogeneous nucleosynthesis could resurrect ~ 1 — the speculations about fun- 
damental physics beyond the standard model do. These well founded speculations include 
Peccei-Quinn (PQ) symmetry, technicolor, supersymmetry, grand unification, and super- 
strings, In the context of the hot big bang model these speculations lead to the prediction 
of various cosmological relics, including particles, topological defects (cosmic strings, do- 
main walls, monopoles, and textures), and the baryon asymmetry of the Universe. The 
discovery — or nondiscovery — of an expected relic provides an important cosmological win- 
dow on fundamental physics beyond the standard model. Since terrestrial experiments are 
hard pressed to probe the physics beyond the standard model, the Heavenly Laboratory 
has become an indispensable testing ground for fundamental physics. 

In the context of the dark matter problem, the implications of theories that go beyond 
the standard model have great significance. Many of these theories predict particle relics 
whose contribution to the present mass density is comparable to the critical density! This 
is no mean feat, and for many of us is a strong hint that the idea of nonbaryonic particle 
relics as the dark matter is on the right track. 

I have organized my discussion of particle dark-matter candidates into six broad cat- 
egories: thermal relics; “skew” relics; axions; nonthermal relics; “significant-other” relics; 
and exotic relics. I have given the axion is own category not just because it is my favorite 
candidate, but also because the story of relic axions is a very rich one and spans three 
categories! 

Thermal Relics 

Because the Universe was in thermal equilibrium at early times essentially all the 
known particles — and perhaps many particles that are yet to be discovered — were present 
in great abundance: When the temperature T was greater than the mass m of a species, a 
number comparable to that of the photons If thermal equilibrium were the whole story, it 
would be a very uninteresting one indeed: At low temperatures the equilibrium abundance 
of a species is exponentially negligible, a factor of order (m/T) 3 ^ 2 exp(— m/T) less than 
that of the photons. 

A massive particle species can only maintain its equilibrium abundance so long as the 
rate for interactions that regulate its abundance is greater than the expansion rate of the 
Universe: T > H , where the expansion rate of the Universe H = l.67gl^ 2 T 2 /mpi ( g » counts 
the total number of degrees of freedom of all relativistic species and mpi = 1.22 x 10 19 GeV). 
The expansio n rate enters because it sets the rate at which the temperature is decreasing, 
H = |T|/T, and therefore the rate at which phase-space distribution functions must change. 

If we specialize to the case of interest for particle dark matter, a stable (or very long 
lived) particle, the reactions that control the abundance are pair production and anni- 
hilation, and their rates are related by detailed balance. The problem now reduces to a 
textbook example! The particle’s number density n is governed by the Boltzmann equa- 
tion, which takes the form 26 

CLTL 

— + 3Hn = -(ctM) ann (n 2 - n| Q ) , (11) 

where (<t|u|)ann is the thermally averaged annihilation cross section times relative velocity 
and ueq is the equilibrium number density. It is more convenient to recast Eq. (11) in 
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terms of the number of particles per comoving volume/ Y = n/s , where s = 2w 2 g m T 3 /45 
is the entropy density, and the dimensionless evolution variable x = m/T : 


dY 

dx 


XS (<7|u|)aNN fv 2 \ 

- H{T = m) [Y ~ Yeq) ’ 


( 12 ) 


where Ieq = 0.278<7 e ff/<7* (for x <C 3) and 0.145(y/y,)x 3 ^ 2 exp(— x) (for i > 3), j is the 
species’ number of internal degrees of freedom, and <jr e ff = 9 (for bosons) or 0.75^ (for 
fermions). Eq. (12) is a particular form of the Ricatti equation that has no closed form 
solutions; it can be solved easily by approximation or numerical integration. I will highlight 
the evolution of a species’ abundance. 

Roughly speaking, the abundance tracks equilibrium until “freeze out,” which occurs 
at temperature Tp, defined by T = H , where T = ^eq^MJann * s ann ihilation rate 
per particle. After that, annihilations cannot keep pace with the decreasing equilibrium 
abundance (“they freeze out”), and thereafter the number of particles per comoving volume 
remains roughly constant, at approximately its equilibrium value at freeze out: Y^ ~ 
Y(Tf)- The mass density contributed by the relic particles today is 


p — mYaoSQ or Clh 2 = O^SY^m/ eV), 


(13) 


where So — 7.1n 7 ~ 2970 cm 3 is the present entropy density. 

Hot and cold relics 

There are two limiting cases: hot relics — species whose annihilations freeze out while 
they are still relativistic (xp < 3); and cold relics — species whose annihilations freeze out 
while they are nonrelativistic (xp > 3). For a hot relic the present abundance is comparable 
to that of the photons, i.e., Y is of order unity. The weak interactions keep ordinary 
neutrinos in thermal equilibrium until a temperature of a few MeV ; thus a neutrino species 
lighter than a few MeV is a hot relic, and 


_ 0.278<7eff 

" g*(T F ) 


~ 3.9 x 10" 2 




92/i 2 eV ' 


(14) 


(There is an intermediate regime, referred to as warm relics; in this case the freeze out 
temperature is sufficiently high so that g,(Tp ) 1 and Y ^ is significantly less than order 

unity. For example, if Tp > 300 GeV, g , is at least 106.75, which is the total number of 
degrees of freedom in the standard model, and for a fermion with two degrees of freedom, 
e.g., a light axino or gravitino, fl = m/910h 2 eV.) 

Freeze out for a cold relic occurs when the species is very nonrelativistic and the 
species’ present abundance is significantly less than that of photons ( < 1). In this 
very interesting case the relic abundance is inversely proportional to the annihilation cross 
section, 


Y n 


4sf/ v / 57 

mmpi^luDANN’ 


(15) 


1 In the absence of appreciable entropy production, the entropy per comoving volume 
S = R 3 s is conserved, implying that s cc R -3 ; thus the number of particles per comoving 
volume N = R 3 n oc n/s. 
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Fig. 3. Upper: Freeze out of a stable, massive particle species. Solid curves indicate 
equilibrium abundance, and broken curves indicate actual abundance (for different values 
of the annihilation cross section). Lower: Contribution of a massive, stable Dirac neutrino 
species to the present mass density as a function of mass. As explained in the text, Clh 2 
increases as m v for m u < MeV; decreases as m~ 2 for m„ > MeV; and increases as m 2 for 
m„ > lOGeV, thereby achieving f lh 2 = 1 for three values of m„. The general behaviour 
of Q/i 2 vs. mass for any stable particle species is similar (e.g., the neutralino). 





where freeze out occurs for i = if — ln[0.04mpim(<r|u|)ANN<7/\/^- In most cases of 
interest freeze occurs at if ~ 20 — 30, corresponding to Tp ~ m/20 — m/30 (in any case 
if only varies logarithmically). 

This is a rather remarkable result: The relic abundance varies inversely with the 
strength of the species’ interactions — implying that the weak shall prevail! Moreover, 
specifying that the species provides the critical density determines the annihilation cross 
section: (ct|v|)ann ~ 10 -37 cm 2 — roughly that of the weak interactions! 

Massive neutrinos 

The simplest example of a cold relic is a “heavy” neutrino (mass greater than a few 
MeV). Provided its mass is less than that of the Z° boson, (<t|v|)ann ~ G 2 F m 2 and 

(^v) -3 ^ *(<&)-■ < i6 > 

That is, the relic abundance of stable neutrino whose mass is a few GeV would provide 
closure density. Since none of the three known neutrino species can be this massive and 
the SLC/LEP results rule out a fourth neutrino (unless it is heavier than about 45 GeV), 
this result, first discussed by Lee and Weinberg, 27 is only an interesting example. 

For a neutrino more massive than about 100 GeV the annihilation cross section begins 
to decrease as m -2 , due to the momentum dependence of the Z° propagator. In this 
regime oc m, and f lh 2 varies as m 2 , increasing to order unity for a mass of order a few 
TeV. 28 

Bringing everything together, the relic mass density of a stable neutrino species in- 
creases as m up to a mass of a few MeV; it then decreases as m 2 up to a mass of order 
100 GeV; and finally it increases as m 2 for larger masses. A stable neutrino species can 
contribute critical density for three values of its mass: C?(100eV); 0(1 GeV); and 0( TeV). 
This behavior is generic for a particle whose annihilations proceed through a massive boson 
(here the Z°). 

Griest and Kamionkowski 29 have generalized this result. Unitarity provides a bound on 
the annihilation cross section of any pointlike species: (<t|u|)ann < Sn/m 2 . This implies a 
lower bound to f \h 2 that increases as m 2 ; requiring that f lh 2 be no larger than one (based 
upon the age of the Universe 30 ) results in an upper bound of 340 TeV to the mass of any 
stable, pointlike species. 

Neutralinos 

A more viable cold relic is the lightest supersymmetric partner or LSP. In supersym- 
metric extensions of the standard model a discrete symmetry, R- parity, is usually imposed 
(to ensure the longevity of the proton); it also guarantees the stability of the LSP. In most 
supersymmetric extensions of the standard model the LSP is the (lightest) neutralino (it 
could in principle be the sneutrino or gluino). The neutralino(s) are the four mass eigen- 
states that are linear combinations of the Bino, Wino, and two Higgsinos. In many models 
discussed early on, especially ones where the LSP was relatively light, the neutralino (by 
which I mean the lightest neutralino) was almost a pure photino state, and thus was 
referred to as the photino. 

The minimal supersymmetric extension of the standard model has a number of param- 
eters that must be specified: jJ. and M, two soft supersymmetry breaking mass parameters 
which are expected to be 100 GeV to few TeV; tan/5 = v 2 /v i, the ratio of the two Higgs 
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vacuum expectation values; the top quark mass; and the scalar quark and scalar lepton 
masses. 31 These parameters determine the composition of the neutralino, its mass, and 
its interactions. The parameter space of supersymmetric models is multidimensional and 
cumbersome to deal with. 

To determine the relic neutralino abundance all one has to do is calculate the cross 
section for neutralino annihilation (the neutralino is a Majorana fermion). For a neutralino 
that is lighter than the W ± boson, the final states Eire fermion-antifermion pairs and light 
Higgs bosons. For the most general neutralino this task has been done by Griest. For 
neutralinos that are heavier than the W ± boson, many additional final states open up: 
W + W~ , Z°Z°, HH, HW, and HZ. This complicated cross section has been calculated 
by Kamionkowski and his collaborators. 33 Let me summarize the salient points. 

• Because the scale of supersymmetry breaking is roughly of order the weak scale, “spart- 
ner” masses are of order the weak scale; since the interactions of the neutralino with 
ordinary matter involve the exchange of spartners, bosons, Higgs bosons, or Z° 
bosons, the neutralino’s interactions are roughly weak in strength. Many of the quali- 
tative features of the relic neutralino abundance are the same as for a neutrino. 

• Over almost the entirety of the parameter space of the minimal supersymmetric exten- 
sion of the standard model the relic neutralino abundance Q. x h 2 is greater than 10 -3 ; 
and in large regions of parameter space D x h 2 is of order unity. This of course traces 
to the fact that the neutralino’s interactions with ordinary matter are roughly weak, 
and makes the neutralino a rather compelling dark matter candidate. 

• Neutralinos can provide the critical density for masses from order 10 GeV to order 
3TeV (depending upon the model parameters). Fixing some of the parameters and 
examining Q, x h 2 as function of m x reveals a similar behavior as for neutrinos: D x h 2 ~ 1 
for a mass in the GeV range and for a mass in the TeV range. 

• Just as with a heavy neutrino, for large neutralino masses the annihilation cross section 
decreases as 1/m 2 ; this results in a maximum neutralino mass that is cosmologically 
acceptable: 3.5 TeV. For m x > 3.5 TeV, Cl x h 2 is greater unity for all models. 

• Finally, the parameter space of models is constrained by unsuccessful accelerator-based 
searches for evidence of supersymmetry. Broadly speaking, the failure to find any 
evidence for supersymmetry has slowly pushed the expected mass of the neutralino 
upward. 34 

Axinos — A Dark Horse LSP 53 

In low-energy supersymmetric models that also incorporate Peccei-Quinn symmetry 
(see Axions below) the axion has a supersymmetric fermionic partner called the ax- 
ino. There are two possibilities for the mass of the axino: (i) of order asm susy ~ 
10 GeV— 100 GeV; or (ii) of order m§ US Y /(f a /N) which is (3(few keV) for f a /N ~ 10 12 GeV 
(msusY ~ 100 GeV — 1 TeV is the scale of supersymmetry breaking). This makes the axino 
a serious candidate for the LSP. In case (i), if the axino is the LSP its relic abundance 
is far too large; even if it isn’t the LSP its decays lead to cosmological havoc, including 
overproduction of the LSP and disruption of primordial nucleosynthesis. Case (i) appears 
to be cosmologically excluded. 

Case (ii) is very intriguing. The axino has a mass in the keV range and is clearly 
the LSP. Such axinos would be brought into thermal equilibrium in the early Universe 
(gluon + gluino — > gluon + axino) and decouple at a temperature of order 10 10 GeV when 
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g*{Tp) ^ 230. Their relic abundance is Yoo = 0.278g e ff/ g+(Tp) ^ 2 x 10 -3 leading to 
ftaxino^ 2 < m ax i no /2keV. That is, for interesting values of msusY and f a /N axinos could 
provide closure density as a warm relic as well as rendering the neutralino impotent. 

How accurately are relic abundances known ? 

Calculating the relic abundance of a species that was once in thermal equilibrium 
has become a routine chore for the particle cosmologist. Because of the importance of 
this calculation, it is prudent to consider the inherent uncertainties. They axe easy to 
identify. 35 Recall that freeze out involves the competition between the expansion rate and 
the annihilation rate. The annihilation rate as a function of temperature is determined by 
the properties of the species — and is thus a given. In calculating the expansion rate we 
have assumed that the Universe was radiation dominated at freeze out; further we assumed 
that there was no entropy production since freeze out, so that remains constant. 

• If there the entropy per comoving volume increased by a factor of 7 after freeze out, 
then the relic abundance Y^ is decreased by the same factor 7. Entropy release could 
occur in a first-order phase transition, or through the out-of-equilibrium decay of a 
massive paxticle species. 

• Additional forms of energy density in the early Universe (e.g., scalar fields, or shear) 
serve to increase the expansion rate at fixed temperature. This in turn leads to an 
earlier freeze out, at a larger abundance. Increasing H(T ) then can increase Yoq. 
While we can be confident that the Universe was radiation dominated by the epoch of 
nucleosynthesis, freeze out for most dark matter candidates occurs earlier, at a time 
when we cannot exclude the possibility that there were additional contributions to the 
energy density. 

Skew Relics 36 

In discussing thermal relics I tacitly assumed that the abundance of the particle and 
its antiparticle were equal. For a Majorana fermion (like the neutralino) this is necessarily 
so; a Dirac fermion (or a scalar species) can carry a conserved (or at least approximately 
conserved) quantum number, and if the net particle number is sufficiently large it will 
determine the relic abundance of the species. Baryon number provides a simple example; 
if there were no net baryon number, baryons and antibaryons would annihilate down to a 
relic abundance n b /s = n b js ~ 10 — 19 , which is significantly smaller than that observed, 
n h /s ~ 77/7 10 10 . As is well appreciated the relic baryon abundance is determined by 

the net baryon number: n^/s = ng/s (the net baryon number density rig = nj — n~ b ). 

The same can occur for any species whose net particle number is conserved, e.g., a 
heavy Dirac neutrino whose net particle number is conserved because of conservation of 
family lepton number. Denote the net particle number per comoving volume by ng/s ( L 
for lepton number). Since the relic abundance cannot be less than the net particle number, 
it follows roughly that: If the net particle number is greater than the would-be freeze out 
abundance, the relic abundance is determined by it, Y 0 0 = ng/s; on the other hand, if the 
net particle number is smaller than the would-be freeze out abundance, the net particle 
number plays no important role and the relic abundance is given by the usual freeze out 
abundance, Yoo ~ Y(x f ). 
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In the case that the relic abundance is determined by the net particle number 



that is, a particle species of mass 35 GeV with a net particle number comparable to the 
baryon asymmetry would contribute the critical density. 

Axions 

Peccei-Quinn (PQ) symmetry with its attendant pseudo-Nambu-Goldstone boson — the 
axion — remains the most attractive and promising solution to the strong-C P problem. 37 
Moreover, the axion arises naturally in supersymmetric and superstring models. One 
might call PQ symmetry and the axion the simplest and most compelling extension to the 
standard model! 

The axion mass and PQ symmetry breaking scale are related by 


Jz 0.62 eV 

777 ~ ~ — — 

a - 1+z (f a /N)~ (f a /N)/W GeV' 


(18) 


where /„ is the PQ symmetry breaking scale, z ~ 0.56 is the ratio of the up to down 
quark masses, f w and m K are pion decay constant and mass, and N is the color anomaly 
of PQ symmetry. At present there is little theoretical guidance as to the key parameter: 
the axion mass, although a variety of astrophysical and cosmological arguments leave open 
only two “windows” for the axion mass: 38 10 -6 eV to 10 -3 eV and 3eV to 8eV (hadronic 
axions only). 

Relic axions arise due to three distinct mechanisms: thermal production 39 — for an 
axion of mass greater than about 10 -4 eV axions thermalize shortly after the QCD tran- 
sition and, today, like neutrinos, should have a relic abundance of order 30 cm -3 ; and 
two coherent processes, the “misalignment” mechanism 40 (see below) and axionic string 
decay 41 — since PQ symmetry breaking involves the spontaneous breakdown of a global 
U(l) symmetry, strings are produced; they decay by radiating (among other things) ax- 
ions. While the thermal population of axions dominates for axion masses greater than 
about 10 -2 eV, there are strong astrophysical constraints in this mass range which pre- 
clude an axion more massive than about 8eV. Thus, thermal axions can contribute at 
most 10% of critical density (more later on thermal axions). 

For axion masses greater than about 10 -2 eV misalignment and axionic string decay 
are the dominant production processes, and sufficient numbers of axions can be produced 
to provide closure density. The importance of axionic string decay is still a matter of 
intense debate. It seems to be agreed that axion production through this mechanism is 
somewhere between being comparable to and about 100 times more important than the 
misalignment mechanism, 41 further that if the Universe inflated either before or during 
PQ symmetry breaking, the number of axions produced by axionic strings is negligible. In 
the “no inflation” case, if axionic string decay is as potent as is claimed by some authors, 
axions provide the critical density for an axion mass of about 10 -3 eV. 

Let me briefly describe the misalignment mechanism. The free energy of the vacuum 
depends upon the axion field because this field modulates the phase of the instanton 
amplitude. At low temperatures the free energy has a maximum value of about Aq CD , is 
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periodic in the “axion angle” 9 = a/(f a /N), and is minimized at a value of 9 = 0. The 
mass of the axion is determined by the curvature of the free energy at 9 = 0 and is given 
approximately by Eq. (18). At high temperatures instanton effects axe strongly suppressed, 
and for T Aqcd the free energy is essentially independent of the axion field. Thus, 
when PQ symmetry breaking occurs (T ~ / a ), no value of the axion angle is singled out 
dynamically, and one expects that the value of the axion angle in different causally distinct 
regions will be randomly distributed between — tt and 7r. Thus the primeval energy density 
associated with the misalignment of the axion field should be of order Aq CD . Around a 
temperature of order Aqcd instanton effects become potent, and the axion mass starts to 
“turn on.” When the axion mass exceeds ZH the axion field will begin to relax toward 
9 = 0. Because it has no efficient way to shed energy, the field is left oscillating. The 
energy density in oscillations of the axion field behaves as nonrelativistic matter during 
the subsequent evolution of the Universe, and may be interpreted in particle language as 
a gas of zero-momentum axions. 

The contribution of these axions to the present mass density of the Universe is esti- 
mated to be 40 

n a h 2 ~ 0.13 x 10 ±OA A^ 0 7 f(9 2 )9 2 (m a /l0~ 5 eV)" 118 . (19) 

where Aqcd = A 2 oo 200 MeV, and 9\ is the initial misalignment angle. The function f(9\) 
accounts for anharmonic effects, and is of order unity (and specifically / — ► 1 for 9\ <C 1). 
The 10 ±0,4 factor is an estimate of theoretical uncertainties — e.g., in the temperature 
dependence of the axion mass. Provided that 9\ ~ 0(1) closure density in axions is 
achieved for a mass somewhere between 10 -6 eV and 10 -4 eV, and for a mass less than 
about 10 -6 eV axions “overdose” the Universe." 1 

The unusual dependence of the axion energy density upon the axion mass is easily 
understood. Regardless of the value of the axion mass, the energy density associated with 
the initial misalignment of the axion field is of order Aq CD ; once the axion field starts to 
oscillate that energy density red shifts as R~ 3 . The axion field begins to oscillate when 
the axion mass m a (T) ~ ZH: For smaller masses the axion oscillations begin later, and 
the energy density trapped in the misalignment of the axion field is diminished less. 

Since the initial misalignment angle 9\ is a random variable, at the time of PQ symme- 
try breaking the value of 9 \ will be different and uncorrelated in different causally distinct 
regions of the Universe. In the absence of inflation, these different regions are very small, 
and today the Universe is comprised of a very large number of regions that each had a 
different value of 9\ . To obtain the average axion energy density, one uses the rms average 
of which is just tt/ 3, in Eq. (19). In this circumstance axions provide closure density 
for a mass in the range of 10 -6 eV to 10 -4 eV. 

If the Universe inflated before or during PQ symmetry breaking the fluctuations in the 
axion field take an entirely different form. While the average of 9\ over many causally- 
separate volumes is still 7t/3, the practical relevance of this fact is nil, because the entire 

m Overdose is not completely accurate; if the Universe is open, the production of axions — 
or any other particle — cannot change the geometry and close it. More precisely, a larger 
value of Slh 2 leads to an earlier epoch of matter-radiation equality and ultimately to a 
more youthful Universe. Requiring that the Universe be at least 10 Gyr old and h > 0.4 
constrains tth 2 < l. 30 
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Fig. 4. Contribution of relic axions to the present mass density as a function of axion mass. 
Subscript “TH” indicates the contribution of thermal relic axions; “MIS” the contribution 
of axions produced by the misalignment process; “S” the contribution of axions produced 
by the decay of axionic strings. Note, in the case that the Universe inflated after or during 
PQ SSB breaking Qs = 0, and I^mis is proportional to the misalignment angle squared, 
whose value is unknown. 




presently observable Universe lies within one causal region where 9 X is constant. A number 
of authors 42 have pointed out that an axion of mass smaller than 10 -6 eV could lead to 

~ 1, provided that 9\ was sufficiently small: 

61 ~ /i(m a /10" 6 eV)°- 59 . (20) 

In this case, then, we would be living in a rare, axion-poor region of the Universe. If the 
Universe did indeed undergo inflation, the fundamental laws of physics do not determine 
9 X . Despite its cosmic import the local value of this parameter is an “historical accident,” 
and can only be determined through direct measurement of fl„/i 2 and m a . n 

Nonthermal Relics 

The axion provides two examples of how a relic can be produced coherently rather 
than thermally: the misalignment mechanism and axionic string decay. For both of these 
processes the number of axions produced is highly superthermal, as is clear since theses 
productions mechanism dominate thermal production for m„ < 10 -2 . 

There axe other examples of nonthermal relics. The most familiar is the superheavy 
magnetic monopole. The monopole is a topologically nontrivial configuration of gauge and 
Higgs fields. Monopoles are produced as topological defects in a symmetry breaking phase 
transition where a semi-simple group Q is broken down to a smellier group 7i that contains 
a 17(1) factor; e.g., SU( 5) — > 517(3) ® 517(2) ® 17(1). Because of the finite size of the 
particle horizon in the standard cosmology, after symmetry breaking the Higgs field can 
only be correlated on distance scales less than H~ l ~ ct — and thus must be uncorrelated 
on larger scales. Because of this fact of order one monopole per horizon volume will be 
produced. Monopole annihilation is ineffective, the monopoles produced should be with 
us today. This production process, which relieves on the fact that the Higgs field cannot 
be correlated on scales larger than the horizon, implies quite generally that order one 
topological defect per horizon volume should arise in a phase transition. It is known as 
the “Kibble mechanism.” 

For the simplest symmetry breaking patterns GUT monopoles are so copiously pro- 
duced by the Kibble mechanism that they overdose the Universe by a factor of about 
10 10 ! Moreover, there axe other stringent astrophysical bounds to their relic abundance. 
Inflation solves the monopole problem by expanding the horizon to a size that is larger 
than our present Hubble volume, and thus predicts less than one monopole in the Universe 
due to the Kibble mechanism. 

More complicated symmetry breaking schemes can reduce the relic monopole abun- 
dance to an acceptable level; and it is possible that significant numbers of monopoles can 
be produced as thermal pairs after inflation. It is very difficult to make a sensible predic- 
tion for the relic abundances of monopoles; however, magnetic monopoles of mass 10 19 GeV 
could provide closure density and have a flux that is consistent with all the astrophysical 
constraints. 44 

n One might then be left with the impression that if the Universe underwent inflation, any 
axion mass can provide closure density provided that 0 X is appropriately small. Additional, 
very important constraints emerge when fluctuations in the axion field that arise during 
inflation are taken into account. 43 
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There are other examples of nonthermal relics, including soliton stars. 45 Soliton stars 
are regions of false vacuum that are stabilized by dynamics rather than topology. (By 
contrast, magnetic monopoles, domain walls, and cosmic string are regions of false vacuum 
that are stable for topological reasons.) For example, imagine a closed region of false 
vacuum associated with a scalar field <f>. Such a region is unstable and should collapse. 
However, if there are particles inside this region whose mass when they are in the false 
vacuum is less than when they are in the true vacuum, they can exert pressure and stabilize 
the region. Whether soliton stars are an interesting dark matter candidate remains to be 
seen. 

“Significant-Other” Relics 

Up to this point I have focused on relics that contribute the critical density. A relic 
from the early Universe can be interesting and significant even if it contributes only a 
fraction of the critical density; e.g., most cosmologists consider baryons (fie ~ 0.1) and 
microwave photons (ft 7 ~ 10 -4 ) to be interesting relics, in spite of their small contributions 
to ft! I will use the term “significant-other” for such relics. 

I will mention two possible significant-other relics: a neutralino and an axion of mass 
3eV to 8eV. While it is possible that the neutralino contributes the critical density, it 
need not be the case. However, in the minimal supersymmetric extension of the standard 
model, the neutralino contributes at least 0.1% of the critical density; thus, if Nature 
exhibits low-energy supersymmetry, the neutralino is at the very least a significant-other 
relic! Moreover, efforts to directly detect relic neutralinos could still be successful even if 
they are only a significant-other relic. Needless to say the implications of their discovery 
for cosmology and particle physics would be almost as profound. 

Axions of mass 3 eV to 8 eV arise as thermal relics and would contribute only about 1% 
of the critical density. Such an abundance is sufficient to permit their detection through 
their decay to two photons. 46 The axion mean lifetime 

r(a — * 27 ) ~ 6.8 x 10 24 £ -2 (m a / eV) -5 sec, (21) 

where £ = [E/N — 2(z+4)/3(.r-|-l)]/0.72 ~ (E/N — 1.95)/0.72 and E is the electromagnetic 
anomaly of PQ symmetry. In the simplest axion models, E/N = 8/3 and £ = 1. 

Relic thermal axions will fall into the various potential wells that develop in the Uni- 
verse as structure formation proceeds. Today they will be found in extended structures such 
as the halos of galaxies and clusters of galaxies, as they cannot dissipate energy and col- 
lapse further. They will decay and produce photons of wavelength A a ~ 24800 A/(m a / eV). 
This radiation will be Doppler-broadened due to the velocities that axions have in these 
objects — for galaxies AA/A ~ v/c ~ 10 -3 and for clusters AA/A ~ v/c ~ 10~ 2 — and for 
distant objects the line will also be red shifted. The most favorable case for their detection 
is to search for the radiation from decaying axions in clusters. The intensity of the axion 
line is approximately 47 

-Cluster ~ 10 — 1 ' £ 2 (m a / 3 eV) 7 erg cm -2 arcsec -2 A -1 s -1 /(l + z c ) 4 , 
where z c is the red shift of the cluster. 

The background against with which this line must compete is the “night sky,” which at 
a ground-based observatory is dominated by the glow of the atmosphere and includes many 
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Fig. 5. Upper: Spectrum of A2218 close to the cluster core (top curve). The spectrum is 
dominated by the “night sky” (atmospheric emission). Spectrum of A2218 near the cluster 
core minus a spectrum at five cluster-core radii from the core (bottom curve). The narrow 
features at 5035A and 5350A are cosmic-ray hits. Lower: “On-off” spectrum for A2256 
with the line expected for a 3.2 eV axion artificially introduced. Note the factor of ten 
change in scale from the upper Figure to the lower Figure. 





strong lines. The baseline intensity of the night sky is 10 -17 erg cm -2 arcsec -2 A -1 s -1 . By 
subtracting “off-cluster” measurements from “on-cluster” measurements one can eliminate 
the night-sky background. This past May, two students, M. Ted Ressell and Matthew 
Bershady, and I used the 2.1m telescope at Kitt Peak to search for axion radiation in 
three clusters using this technique. The spectra we took span 3600 A to 8600 A with 
10 A resolution. Our “on-off” subtractions allowed us to search for such a line with a 
sensitivity of less them 3% of the night sky for the mass range from 3.1 eV to 7.9 eV. 
Unfortunately, our search proved unsuccessful, and we have closed this mass window. 48 
(The lower mass limit to this window, 3eV, derives from the SN 1987A limit. Obviously 
there are uncertainties inherent in this limit, and perhaps an axion of mass 2eV to 3eV 
is still permitted. Atmospheric emission — OH bands — preclude a ground-based search for 
such an axion; however, one could search the mass range of 2 eV to 3 eV using the Hubble 
Space Telescope! A proposal is in the works.) 

Exotic Relics 

Thus far I have focused on particle relics that today would behave like ordinary nonrel- 
ativistic matter. There are more exotic possibilities. Since the amount of matter associated 
with bright galaxies seems to contribute only 20% of critical, and a strong theoretical prej- 
udice for ft — 1 exists, several relics have been suggested that today would contribute an 
almost uniform energy density of 80% of the critical density. A uniform contribution to the 
mass density would not show up in the dynamical measurements, thereby solving the “ft 
problem.” The exotic candidates include a relic cosmological constant, 49 very light cosmic 
strings that are either fast moving or exist in a tangled network, 50 or relativistic particles 
produced by the recent decays of a massive relic. 51 Whether or not we have to resort to 
such exotics to savage our strong prejudice remains to be seen. 

IV. Implications for Structure Formation in the Universe 

According to the standard cosmology, structure formation proceeds via the Jeans (or 
gravitational) instability: Small primeval density perturbations begin to grow once the 
Universe becomes matter dominated, and then develop into the structure that we observe 
today. The structure-formation problem is essentially an initial data problem: Specify the 
primeval density perturbations and the quantity and composition of matter, and let it go! 

We now have well motivated suggestions for both pieces of initial data. 52 For the density 
perturbations, there are several choices: inflation-produced, constant-curvature (Harrison- 
Zel’dovich) perturbations; inflation-produced, isocurvature perturbations; and topological 
relics, such as cosmic strings or texture, as the seed perturbations. For the matter content, 
there are the following suggestions: ft = 1 , ftg ~ 0.1, and ftx ~ 0.9, where generically X is 
hot dark matter (a light neutrino species), cold dark matter (axions, neutralinos, magnetic 
monopoles, ...), or perhaps warm dark matter (a keV mass particle, such as an axino, 53 
righthanded neutrino, or gravitino). 0 

0 In the context of structure formation, hot, warm, and cold refer to the velocity disper- 
sion of the relic particles around the time of matter-radiation equality; hot corresponds to 
relativistic and cold to very nonrelativistic. For thermal relics this matches the previous 
nomenclature; in general nonthermal relics have very small velocity dispersions and behave 
like cold dark matter. 
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The suggestion that weakly interacting relic particles comprise the bulk of the mass 
density of the Universe and contribute fi ~ 1 has been a particularly important one, and 
virtually all scenarios of structure formation now include nonbaryonic dark matter. For 
good reason; in a “particle dark-matter” Universe density perturbations can begin growing 
as soon as the Universe becomes matter dominated, while in a baryon-dominated Universe 
density perturbations cannot begin to grow until decoupling. Further, linear perturbations 
in a low-fi model cease growing at a red shift z ~ fi -1 . Thus, in a low-fi model larger 
amplitude perturbations are required. Low-fi models with curvature-perturbations conflict 
with the observed isotropy of the CMBR if fi < 0.3. 

Two “stories” of structure formation have been studied in some detail: hot dark matter 
and cold dark matter (both with inflation-produced, constant-curvature perturbations). 
Hot dark matter seems to be ruled out, as galaxies form too late. 54 Cold dark matter 
is the most successful paradigm for structure formation yet proposed. 55 Other scenarios 
involving cosmic strings and texture are presently less well developed. In any case, the 
“hints from the early Universe” as to the initial data for structure formation have served 
well to bring this problem into sharper focus. Next, I will digress briefly to discuss my 
candidate for the “best-fit model” of the Universe. 


The Best-fit Universe 

Cold dark matter does a remarkably good job of describing the Universe on scales less 
than about 20 h~ l Mpc. However, it appears to have a number of shortcomings: deficient 
large-scale structure, deficient galaxy counts, the age problems, and the fi problem. No 
one of these problems is sufficiently troublesome to falsify the cold dark matter paradigm — 
yet — but taken together they are worrisome. As we shall see, the addition of a cosmological 
constant simultaneously addresses all of these problems. 

As a reference point, the conventional cold dark matter scenario is: a flat Universe 
whose composition is fie ~ 0.1 <C ficDM ~ 0.9, with h ~ 0.5 (to have a sufficiently 
old Universe) and inflation-produced Harrison-Zel’dovich curvature perturbations whose 
spectrum after the epoch of matter-radiation equality is 55 


h | 2 


Ak 

(l + jSJfe+w* 1 - 5 +7ifc 2 ) 2 ' 


( 22 ) 


Here 6 k is the amplitude of the Fourier component of comoving wavenumber k (= 27r/A), 
A is an overall normalization constant, 0 = 1.7(fi/t 2 ) -1 Mpc, u = 9.0(fi/i 2 ) -1 5 Mpc 1 ' , 
and 7 = 1.0(fi h 2 )~ 2 Mpc 2 . 

The basic idea of the best-fit model is simple; retain the flatness, but add a cosmological 
constant. 49,56 The model I discuss here is: (i) Hubble constant of around 70kms -1 Mpc -1 
( h = 0.7) — a nice compromise value; (ii) fig ~ 0.03 — near the central value implied by nu- 
cleosynthesis; (iii) ficDM ~ 0.17 — sufficiently greater than the baryonic component so that 
the mass density is dominated by that of the cold dark matter; (iv) fi A — cosmological 
constant corresponding to an energy density p\ = fiaPcRiT — 3.2 x 10 -47 GeV 4 = 
(2.4 x 10 -3 eV) 4 . I am not wed to these particular values and I simply use this set for 
definiteness. (If the ratio of the mass densities of CDM and baryons is somewhat smaller, 
then the decoupling of matter and radiation can have an effect on the spectrum of density 
perturbations, which is to boost power on large scales. 57 If the “best-fit model” is still 
deficient in large-scale power, this effect could improve the situation.) 
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For this model the total matter contribution Anr = 0.2, and today the vacuum energy 
density dominates the matter energy density by a factor of four. In general the ratio 
PHr/pa = 0.25(1 + z) 3 . At red shifts greater than about z A ~ 0.59 the matter energy 
density dominates, and the model behaves just a flat, CDM model. To determine when this 
model becomes matter dominated one simply sets f 1 h 2 = Anr h 2 ^ 0.098: Teq = 0.54 eV; 
< E q ~ 4.5 x 10 12 sec; and z E q — 2300. Once the radiation energy density is negligible 
( z <C zeq), the scale factor evolves as 

m = (^) ' sinh !/3 ( 3yfthHot/2 ) , (23) 


where the value of the scale factor today is taken to be one. 

The A problem 

A cosmological constant behaves just like a uniform mass density (with equation of 
state p = —p). As such, it would not affect determinations of A based upon dynamics 
(galactic halos and cluster virial masses). These measurements of the masses of tightly 
bound systems are insensitive to the contribution of a uniform background energy density 
because the average density in these objects is much greater than the average density of 
the Universe. Likewise, determinations of A based upon the peculiar velocities induced by 
the clumpy matter distribution would only reveal the clumpy, matter component. Thus, 
all current dynamical determinations that indicate A ~ 0.1 —0.3, would be consistent with 
a flat Universe (A = 1) with Anr = 0.2. 

The age problems 

As is well appreciated the addition of a cosmological constant increases the age of a 
flat Universe. The age of a A model is 


*(*)- 


2V 

3v/Aa 


sinh 1 


\/A a /A nr /(1 + z) 3/2 ; 


(24a) 
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2ff 0 ~ 1 r i + vgr 

3-v/Aa v^Anr . 


(246) 


The present age of a A-model is always greater than 2 H 0 1 / 3 and for A a = 0.8, t 0 = 
Uflo' 1 ~ 15.5 Gyr, an age which is comfortably consistent with the age as determined 
from the radioactive elements, from the oldest globular clusters, and from white dwarf 
cooling (e.g., see Ch. 1 of Ref. 52 and references therein). Moreover, a A model is older 
than its matter-dominated counterpart at any given epoch, so that objec ts at a given red 
shift have had more time to evolve. For z >• z A , t(z) — *• 2^J" 1 /3\/Anr(1 + z) 3 ? 2 , which 
is a factor of A^/ 2 older than a flat, matter-dominated model; at these early epochs the 
“best-fit model” is a factor of 1.6 older than the conventional CDM model. 


Large-scale structure 

The spectrum of density perturbations at matter- radiation equality, ( SM/M ) cc 
Jt 3 / 2 |5fc|, decreases monotonically with A and its wavelength scale is determined by the 
value of A h 2 . The spectrum “shifts” to larger length scales as A/i 2 is decreased. Supposing 
that the spectrum is normalized on the scale A = 8/i -1 Mpc (a common normalization 
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Fig. 6. The spectra of rms mass fluctuations for the best-motivated (conventional CDM) 
and the best-fit models. The rms mass fluctuation is computed using the “top-hat” win- 
dow function (sharp sphere of radius R) and normalized to 6M/M = 1 on the scale 
R ~ Sh~ 1 Mpc. The spectrum for the best-fit model is shifted to the right (relative to 
conventional CDM) because is smaller; for this reason it has more power on large 

scales. 


is: SM/M — 1 for A ~ Sh~ l Mpc), decreasing f ~lh 2 increases the power on all scales 
greater than the normalization scale. Put another way, the ratio of the characteristic 
scale in the spectrum, Aeq = 13(f)h 2 ) -1 Mpc, to the scale of nonlinearity in the Universe, 
Ajjl — 8h -1 Mpc, is Aeq/Anl — 1.6/fiA; in the “best-fit model” this ratio is a factor of 
3.5 greater than in a model with Cl = 1 and h = 0.5 (conventional cold dark matter, or the 
“most well motivated model”), implying more power on large scales. Needless to say, this 
can only help with the problem of deficient large-scale structure. 

To be specific, if the spectrum of perturbations is normalized by (6M/M)x=ah- 1 Mpc = 
l, p I find that: A = 4.4 x 10 6 Mpc 4 for Cl = 1 and h = 0.5 (conventional CDM) and 
A = 2.5 x 10 7 Mpc 4 for = 0-2 and h = 0.7 (“best-fit model”). On large scales 
(A Aeq) 6M/M oc \J~Aj A 2 ; it follows that 8M/M for the “best-fit model” is a factor of 
4.7 bigger on large scales. 

Growth of density perturbations 

Subhorizon-sized, linear density perturbations grow as the scale factor during the 
matter-dominated regime (z < zeq — 23000f2/i 2 ), and remain roughly constant in ampli- 
tude when the Universe is radiation dominated, curvature dominated (z < zcurv — Cl -1 — 
2; -^curv — 3 for ft = 0.2), or vacuum-energy dominated (z A ~ [ftj^ 1 — l] 1 / 3 — 1 ~ 0.59). 
For a nonflat, Cl = 0.2 model the reduction in the growth of perturbations relative to a flat 
model is very significant: about a factor of 20. By contrast, in flat-A models perturbations 
grow almost unhindered until the present (see Refs. 49 and 58). In the “best-fit model” 
the growth factor is only a factor of 0.8 less than zeq, or about 1800. For comparison, in 
the conventional CDM model the growth factor zeq — 5800, only about a factor of three 
more growth. 

Microwave anisotropies 

For conventional CDM the predicted CMBR temperature anisotropies are about a 
factor of three or so below the current level of observed isotropy (depending upon the 
angular scale and biasing factor b). 59 One might worry that because the “best-fit model” 
has more power on large scales and the growth factor for perturbations is smaller the 
predicted CMBR anisotropies might violate current bounds. That is not the case. The 
reason involves the angular size on the sky 9 of a given scale A at epoch z: 


9( A,z) = A/r(z); 

f to JL _ 2g 0~ x r alnh ~ l [\/«a/«Nr| du 

Jt(z) R(t) 3Q^ 6 Qj/^ J sinh -1 [v^nA/fl-t-*)^™] sinh 2 / 3 U 1 

where r(z) is the coordinate distance to an object at red shift z. In a flat, matter- 
dominated model r(z) = 2 [l — l/\/l + zj — ► 2 Hq 1 for z ^ 1, and #(A,z 1) — 

34.4" (A//i -1 Mpc). For the “best-fit model” r(z 1) ~ 3.9Hq 1 and 0(\,z >■ 1) ~ 
17.7" (A/A- 1 Mpc). 

In a flat A-model the horizon is further away and a given length scale has a smaller 
angular size. Since the temperature fluctuations on a given angular scale are related to 

p I have used the “top hat” window function [IU(r) = 1 for r < r 0 and = 0 for r > tq] 
to define M, so that ( 8M/M ) 2 = ( 9 / 27 T 2 ) / 0 °° A: 2 |djt| 2 [sin(A:ro)/A: 3 ro - cos(fcr 0 )/fc 2 r 2 ] 2 dk, 
where ro = S h~ x Mpc. 


(25a) 
(25 b) 
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the density perturbations on the length scale that subtends that angle at decoupling, 
in the “best-fit model” temperature fluctuations on a given angular scale are related to 
density perturbations on a larger scale A. While the “best-fit model” has more power on 
a fixed (large) length scale, a fixed angle 9 corresponds to a larger length scale, where the 
amplitude of perturbations is smaller because SM/M decreases with A. 

Consider the temperature fluctuations on large-angular scales ( 9 1°); they arise due 

to the Sachs- Wolfe effect and ( ST/T)$ ~ {8p/ p) hor/2 on the scale A(0) when that scale 
crossed inside the horizon. For the Harrison-Zel’dovich spectrum the horizon- crossing 
amplitude is constant, so that 8T /T is independent of angular scale (for 9 >■ 1°). The 
CMBR quadrupole anisotropy is related to the amplitude of the perturbation that is just 
now crossing inside the horizon: A H or ~ 2 H^ 1 ~ 12000 Mpc (conventional CDM) and 
Ahor ~ 3.9Hq 1 ~ 16700 Mpc (“best-fit model”). Evaluating the normalized spectra on 
these scales it follows that the large-angle temperature fluctuations in the “best-fit model” 
are only a factor of 1.2 larger than for conventional CDM, in spite of the fact that the 
“best-fit model” has significantly more power on large scales. 

The amplitude of the temperature fluctuations on small angular scales ( 9 <C 1°) is 
proportional to the amplitude of the density perturbations at the time of decoupling 
(^dec ~ 1000), on the scale A(0). In the “best-fit model” perturbations have grow by a 
factor of about 0.8zdec since decoupling, while those in the “most well motivated model” 
have grown by a factor of zeq- On the other hand the length scale corresponding to the 
angular scale 9 is larger for the “best-fit model.” The net result is that the temperature 
fluctuations on an angular scale of 1° are also only about a factor of 1.2 larger. 

Galaxy counts 

Because the coordinate distance to an object of given red shift is greater in a flat 
A model, there is greater volume per red shift interval per solid angle, which increases 
the number of galaxies in dzdu. To see roughly how this goes, consider the deceleration 
parameter 

9 o = fi(l + 3p/p)/2 = (l-3fl A )/2~-1.2, (26) 

where is the total energy density p divided by the critical energy density and p is the 
total pressure. From Eq. (2) one can see that the galaxy-number count is significantly 
increased by the addition of a cosmological constant, dNoAh/dz = z 2 yiq a l [1 — 3z + • • •] 
compared to z 2 uq al[ 1 + 0.4z + •••]. 

Large-scale motions 

The rms peculiar velocity of a volume defined by the “window function” W(r), averaged 
over all such volumes in the Universe, is 

1 f°° 

{v) = 2 7*1 * 2 KI 2 |W"WI 2 <tt. (27) 

Using a gaussian window function [W ro (r) = exp(-r 2 /2r 2 )] and normalizing the spectrum 
as above, the rms peculiar velocity expected on the scale r 0 = 50 h~ l Mpc is 60 

U50 ^ 83/i -0 ' 9 kms~ l ~ 160 kms -1 (fl = 1, h = 0.5); 
u 5 o ^ 83O- 9 - 33 /i- 0 - 9 kms _1 ~ 200 kms -1 (fi NR , h = 0.7). 
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While the rms peculiar velocity on the scale of 50Mpc is still fax short of 700 km s 1 , it is 
larger, owing to fact that there is more power on large scales. 9 

Motivation 

As its name suggests, it is a model motivated by observations and not aesthetics: 
Conventional cold dark matter is clearly better motivated. In this regard one should keep in 
mind the words of Francis Crick: “Any theory that agrees with all the data at a given time 
must be wrong!” While the conventional CDM model has one question to answer — why 
the ratio of the baryon density to that of cold dark matter is of order unity (see below) — in 
the “best-fit model” one must also address “why now?” — why is the cosmological constant 
just now becoming dynamical important? (This problem is similar to the flatness problem, 
where the question is, why is the curvature radius just now becoming comparable to the 
Hubble radius?) Moreover, there is the issue of the cosmological constant itself: At present 
there is every reason to expect a cosmological constant p\ = A/8 7rG ~ m pi 4 * that is some 
122 orders of magnitude larger than observations permit*" (Supersymmetry might be able 
to help in this regard, reducing the estimate to p\ ~ Gp 2 * , which is only 56 orders of 
magnitude too large!) The strongest statement that one can make in defense of a relic 
cosmological constant of the desired size is that no good argument exists for excluding it! 


V. A New Dimensionless Cosmic Ratio 63 

Dimensionless numbers play a crucial role in physics and in cosmology, and attempts 
to understand their origin often lead to important insights. There are a number of dimen- 
sionless ratios in cosmology: the baryon-to-photon ratio, the fractional abundances of the 
light elements, the amplitude of the primeval density perturbations, and the ratio of the 
neutrino and photon temperatures. If there is a significant amount of nonbaxyonic matter 
in the Universe, we have a new dimensionless ratio to understand 

(28) 



In particular we can ask why r is order unity, and not say 10 20 or 10 2 °? 

We can try to express r in terms of fundamental quantities. To begin, write 

(28') 

One of the great successes of particle cosmology is the dynamical explanation of the baryon 
asymmetry, or baryogenesis. 64 While the specific details of baryogenesis are still lacking, 
generally one expects that ng/s ~ e/g», where g * ~ 100 — 1000 counts the number of 
degrees of freedom at the epoch of baryogenesis (10 14 GeV?) and e ~ 10 -8 — 10 -7 is a 
measure of the C, CP violation in the baryon number violating sector and — on general 

9 The comparison of theoretical expectations to the peculiar-velocity data is far more 

complicated than just computing ( v 2 ) for a gaussian window function. 61 The point I wish 

to make here is that adding a cosmological constant increases peculiar velocities. 

r There is one interesting explanation of why the cosmological constant is “probably” 

zero: Coleman and others 62 have argued that due to wormhole effects the wavefunction of 

the Universe is very sharply peaked at zero cosmological constant. 
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grounds — is order (ot/ir) N ( a = <7 2 /47t, g is a Higgs coupling). The quantity nx/s is the 
relic abundance of X particles per comoving volume, a quantity that can be calculated as 
we seen above. Now consider the implications of r ~ 0.1 for the various relics previously 
discussed. 

Heavy neutrino/neutradino 

For a thermal relic like a heavy neutrino or neutralino whose interactions axe weak, 
n x /s l/m 3 x mp\G 2 F (see Thermal relics above). The condition that r be of order unity 
implies 

Gn 

Gp y g*mpi 3 

and thus is related to the fact that the weak scale is much smaller than the Planck scale. 
Axion 

The relic abundance of axions can be expressed as n a /s ~ /a /Aqcd^pi- The condition 
that r be of order unity implies 



A 

mpi 


f— V 

\Aqcd/ 


~ e. 


(30) 


and thus is related to the fact that the PQ symmetry breaking scale is somewhat less than 
the Planck scale. 

Light neutrino 

The relic abundance of a light neutrino species, n„/s, is of order unity. If we assume 
that light neutrino masses arise through the see-saw mechanism, then m„ ~ m^/M, where 
nif is a typical fermion mass and M is the laxge energy that characterizes lepton number 
violation. The condition that r be of order unity implies 


m 2 f/m B ~ eM/g m , (31) 

and thus is related to the fact that fermion masses are much smaller than the scale of 
lepton number violation. 

Skew relic 

Consider a skew relic whose net particle number per comoving volume is comparable 
to that of baryon number (perhaps its net particle number was produced at the same time 
as the baryon number, e.g., a heavy neutrino). In this case the fact that r is of order unity 
is related to the fact that the mass of the skew relic is comparable to that of a nucleon. 

In a sense, all of these relations only tell us what we already knew and put in. However, 
this exercise does illustrate the fact that r can be related to fundamental quantities in 
physics, and raises the hope that this very important dimensionless cosmological ratio 
may some day have a more fundamental explanation. Apparently, that explanation will 
have to wait until we have a better understanding of the various energy scales that arise 
in particle physics. 


VI. Summary 

What do we know about the quantity and composition of the matter in the Universe? 
Most of the matter in the Universe is dark, with luminous matter contributing less than 
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1% of the critical density. The best estimates of the amount of matter associated with 
bright galaxies is Oabg — 0.1 —0.3; however, there are some observations that suggest that 
Q might be larger, perhaps even equal to one. Based upon primordial nucleosynthesis, we 
can be confident that baryons contribute between 1.1% and 12% of the critical density — 
more than that of luminous matter, but far less than the critical density. It is no means 
impossible that baryons account for the entire mass density of the Universe. 

While there may already be evidence for nonbaryonic matter — if fl is indeed 0.2 — if 
our strong theoretical prejudice for = 1 is correct, nonbaryonic matter must account 
for the bulk of the mass density in the Universe. In any case, it is certainly a hypothesis 
worthy of careful consideration. 

Theories of fundamental physics that go beyond the standard model have profound 
implications for the earliest moments of the Universe; indeed, many of us believe that the 
“blueprint” for the Universe traces to events that took place during that epoch. Theories 
that unify the particles and interactions predict the existence of new, stable particles (or 
additional properties for known particles, e.g., neutrino masses), and remarkably enough, 
the relic abundances calculated for a number of these new particles is comparable to that 
required to close the Universe. For many, this is what makes the particle dark-matter 
hypothesis so compelling. Needless to say, the discovery of such a relic would not only 
solve a cosmological puzzle, but would also shed light on the theory that unifies the forces 
and particles. 

By now there is a virtual zoo of particle dark-matter candidates. 65 However, three 
candidates are particularly well motivated and attractive. They are an axion of mass 
10 -6 eV to 10 -4 eV, a neutrino of mass 92h 2 eV, and a neutralino of mass 10 GeV to 
3TeV. Peccei-Quinn symmetry and its axion resolve a nagging and serious difficulty of 
the standard model: the strong-CP problem. The neutralino is a very robust prediction of 
theories that incorporate low-energy supersymmetry. Low-energy supersymmetry provides 
some understanding of the hierarchy problem (the large disparity between the weak scale 
and the Planck scale), and is further motivated by superstring theories. Neutrinos actually 
exist — and come in three flavors! — and in many extensions of the standard model small 
neutrino masses are predicted. Moreover, the first results of the SAGE experiment, 66 
together with the results of the Homestake and Kamiokande II solar neutrino experiments, 
suggest that nonadiabatic MSW neutrino oscillations may be the solution to the solar 
neutrino problem. 67 If this is so, it implies a mass for the fiorr neutrino in the range 
10 -4 eV to 10 -2 eV. Speculating (upon supposition to be sure) that this is the mass for 
the fj. neutrino, a simple see-saw scaling estimate for the r neutrino mass might just put 
it in the cosmologically interesting range. While cold dark matter provides a far more 
promising paradigm for structure formation than does hot dark matter, I am certain that 
cosmology could learn how to live with a neutrino-dominated Universe. 3 4 

Particle dark matter is an attractive and compelling hypothesis, and the next step is 
to test it. A variety of experiments are underway, and more are planned. 68 The experi- 
mental efforts encompass a diversity of approaches, involving conventional laboratory and 
accelerator experiments, large- underground detectors, and experiments built expressly to 

3 Dennis Sciama has recently touted the multitude of astrophysical virtues of an unstable 
neutrino, and has gone so far as to precisely “predict” its mass, 28 eV to 30 eV, and lifetime, 

r = 2 ± 1 x 10 21 sec. 67a 
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detect the dark matter particles in our local neighborhood. The search for evidence of 
supersymmetry is going on at accelerator laboratories all over the world. Indirect evidence 
for the existence of particle dark matter in our own halo could come from the annihila- 
tion products of particle dark matter in the halo or from particle dark matter that has 
accumulated in the sun or earth. The GALLEX and SAGE experiments may well provide 
information about neutrino masses, and a nearby supernova or a long-baseline neutrino 
oscillation experiment could provide definite evidence for neutrino masses. The MACRO 
experiment in the Gran Sasso Laboratory is operating and can search for both relic mag- 
netic monopoles and high-energy neutrinos from particle dark-matter annihilations in the 
sun or earth. First-generation Sikivie-type detectors to search for cosmic axions have been 
built and successfully operated; 69 a second generation detector with sufficient sensitiv- 
ity to detect halo axions in the our neighborhood has been proposed. 70 Low- background, 
cryogenic detectors designed to detect the keV energies deposited by halo neutralinos that 
elastically scatter within the detector are under development in laboratories all over the 
world, and low-background ionization detectors have already been used to search for heavy 
neutrinos and cosmions. 71 

The answer to the simple question — What is the Universe made of? — may well be 
answered soon. If the bulk of the matter in the Universe is nonbaryonic, this discovery 
will rank as one of the most important of the century, and will have profound implications 
for both cosmology and particle physics. 
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